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We report on the event structure and double helicity asymmetry (AlL) of jet production in lon- 
gitudinally polarized p+p collisions at -^=200 GeV. Photons and charged particles were measured 
by the PHENIX experiment at midrapidity |rj| < 0.35 with the requirement of a high-momentum 
(> 2 GeV/c) photon in the event. Event structure, such as multiplicity, pr density and thrust 
in the PHENIX acceptance, were measured and compared with the pythia event generator and 
the GEANT detector simulation. The shape of jets and the underlying event were well reproduced 
at this collision energy. For the measurement of jet All, photons and charged particles were 
clustered with a seed-cone algorithm to obtain the cluster pr sum (pr co ). The effect of detector 
response and the underlying events on p!j? c ° was evaluated with the simulation. The production 
rate of reconstructed jets is satisfactorily reproduced with the NLO pQCD jet production cross 
section. For 4 < p^ c ° < 12 GeV/c with an average beam polarization of (P) = 49% we mea- 
sured All = -0.0014 ± 0.0037 stat at the lowest p^ co bin (4-5 GeV/c) and -0.0181 ± 0.0282 stat at 
the highest p^ c ° bin (10-12 GeV/c) with a beam polarization scale error of 9.4% and a pt scale 
error of 10%. Jets in the measured range arise primarily from hard-scattered gluons with 

momentum fraction 0.02 < x < 0.3 according to pythia. The measured All is compared with 
predictions that assume various AG(x) distributions based on the GRSV parameterization. The 
present result imposes the limit —1.1 < f Q2 dxAG(x, fi 2 = lGeV 2 ) < 0.4 at 95% confidence level 
or J ° 3 2 dxAG(x,fi 2 = lGeV 2 ) < 0.5 at 99% confidence level. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

The motivation of this measurement is to understand 
the spin structure of the proton, particularly the contri- 
bution of the gluon spin (AG) to the proton spin. The 
proton spin can be represented as 



1 

2 proto 



AG + L a + L 



(1) 



where AG is the gluon spin, i.e. the integral of the po- 
larized gluon distribution function, AG = f Q dxAG(x), 
5^ Aq is the quark spin, and L q and L g are the orbital 
angular momenta of quarks and gluons in the proton. It 
was found by the EMC experiment at CERN in 1987 that 
the quark spin contribution to the proton spin is only 
(12 ± 9 ± 14)% After the EMC experiment many 

deep inelastic scattering (DIS) experiments have been 
carried out to measure ^ Aq more precisely. The recent 
analysis by the HERMES experiment [3j reported that 
£ Aq = 0.330 ±0.0n(theo.)±0.025(exp.)±0.028(evol.) 
at a hard-scattering scale /i 2 ~ 5 GeV 2 , which is only 
about 30% of the proton spin. Consequently, the major- 
ity of the proton spin should be carried by the remaining 
components. 

Jet production from longitudinally-polarized p+p colli- 
sions is suited for the measurement of AG because gluon- 
involved scatterings, such as q+g — > q+g or g+g — > g+g, 
dominate the cross section. The double helicity asymme- 
try 



All = 



(2) 



'Deceased 

^PHENIX Spokesperson: jacak@skipper.physics.sunysb.edu 



is the asymmetry in cross section between two beam he- 
licity states. In the All measurement, many systematic 
errors cancel out so that high precision can be achieved. 

Another motivation of this measurement is to study 
the event structure of p + p collisions. A high-energy 
p + p collision produces not only hard scattered par- 
tons but also many particles that originate from soft 
interactions which we call the 'underlying event'. The 
pythia event generator phenomenologically models the 
underlying event on the Multi-Parton Interaction (MPI) 
scheme 4], and can reproduce the event structure oip+p 
collisions measured by the CDF experiment at \fs = 1.8 
TeV We present measurements of event structure at 
lower collision energy, *Js = 200 GeV, and compare them 
with those simulated by pythia in order to examine the 
validity of the PYTHIA MPI scheme. One of the goals of 
the PHENIX experiment at the Relativistic Heavy Ion 
Collider (RHIC) is the determination of AG. PHENIX 
has published results on single particle production; the 
All of 7T° production was reported in [a, 0|- This pa- 
per reports a measurement of jet production. For AG, 
it is valuable to determine the parton kinematics follow- 
ing the collision in order to better control the x range. 
In this work we reconstruct jets, observing a larger frac- 
tion of the partons momentum. This allows improved 
reconstruction of the original parton kinematics and bet- 
ter statistical accuracy for higher x gluons. Since 7r°'s in 
p+p collisions are produced via jet fragmentation, the 
measurements of jet and ir° with same data set have a 
statistical overlap. The size of the overlap was estimated 
to be 40-60% depending on the jet px- Even in such 
overlapped events, measured px of jets does not correlate 
with that of 7r°s, and thus the two measurements have an 
independent sensitivity on x. The fraction of q + g sub- 
process is larger than q + q and g + g subprocesses in the 
present jet measurement, making it sensitive to the sign 
of AG. The STAR experiment at RHIC is also measuring 
inclusive jets to determine AG 0]. These measurements 
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have different types of systematic uncertainties and thus 
one can provide a systematic check for the other. 

The remainder of this paper is organized as follows. 
In Section II, the parts of the PHENIX detector that is 
relevant to the jet measurement are described. In Sec- 
tion III, analysis methods such as particle clustering and 
simulation studies are discussed. In Section IV, results 
on event structure, jet production rate and beam-helicity 
asymmetries are shown. 



II. EXPERIMENTAL SETUP 

The PHENIX detector Q can be grouped into three 
parts; the Inner Detectors, the Central Arms and the 
Muon Arms. The schematic drawing of the PHENIX 
detector is shown in Fig. [TJ In this measurement, the 
Central Arms were used to detect photons and charged 
particles in jets, and the Inner Detectors to obtain the 
collision vertex and beam luminosity. 



PHENIX Detector 




West 



Beam View 



East 




South Side View North 

FIG. 1: (color online) PHENIX detector. 

A. Inner detectors 



The Inner Detectors include the Beam-Beam Counters 
(BBC) and the Zero-Degree Calorimeters (ZDC). 



The BBC is composed of two identical sets of counters 
placed at both the north and south sides of the colli- 
sion point with a 144 cm distance [loj . Each counter 
is composed of 64 sets of PMT plus a 3-cm quartz 
Cerenkov radiator. The BBC covers a pseudorapidity of 
3.0 < I77I < 3.9 over the full azimuth. The BBC measures 
the number of charged particles in forward and backward 
regions to determine the collision time, collision z-vertex, 
and beam luminosity. The timing and z-vertex resolution 
in p+p collisions are about 100 ps and 2 cm, respectively. 

The ZDC is comprised of two sets of hadronic calorime- 
ters placed at the north and south sides of the collision 
point with a 18 m distance [ll|. It covers a 10 cm x 
10 cm area perpendicular to the beam direction, which 
corresponds to 2.8 mrad when viewed from the collision 
point. It consists of alternating layers of tungsten ab- 
sorbers and sampling fibers, and is 150 radiation lengths 
and 5.1 interaction lengths in depth. It measures neu- 
trons in forward and backward regions and is used as a 
local polarimctcr which assures that the beam polariza- 
tion is correctly longitudinal or transverse at the interac- 
tion region by observing the left-right asymmetry in the 
p + p —> neutron + X scattering cross section (l2l [j~3l]. 



B. Central Arms 

The Central Arms consist of a tracking system and 
an electromagnetic calorimeter (EMCal). Pad cham- 
bers (PC) and drift chambers (DC) were used to detect 
charged particles in jets, and the EMCal was used to 
detect photons in jets. 

The EMCal system 14{ is located at a distance of 
5 m from the interaction point. The system consists 



of four sectors in each of the East and West Arms, 
and each sector has a size of 2x4 m 2 . The system is 
composed of two types of calorimeter, lead scintillator 
(PbSc) and lead glass (PbGl). One PbSc module has a 
size of 5.5x5.5x37.5 cm 3 corresponding to 18.0 radiation 
lengths. One PbGl module has a size of 4.0x4.0x40.0 
cm 3 corresponding to 14.4 radiation lengths. The energy 
resolution is ~ 7% at E = 1 GeV. 

The DC system [l5| is located in the region from 2 to 
2.4 m from the interaction point to measure the position 
and momentum of charged particles. The DC system 
consists of one frame in each of the East and West Arms. 
Each chamber has a size of 2.5 mx90° in z-<f> direction 
with cylindrical shape, and is composed of 80 sense planes 
with a 2-2.5 cm drift space in the <f> direction. Each sense 
plane has 24 wires, which precisely measure r-<p position, 
and 16 tilted wires, which measure z position. 

The PC system [15| is composed of multi-wire propor- 
tional chambers in three separate layers, which are called 
PCI, PC2 and PC3, of the Central Arms tracking sys- 
tem. The PCI is located behind the DC and is used for 
pattern recognition together with the DC by providing 
the z coordinate. The PCI consists of a single plane of 
anode and field wires lying in a gas volume between two 
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cathode planes. One cathode is segmented into pixels 
with a size of ~ 8.5 x 8.5 mm 2 , and signals from the 
pixels are read out. 

Charged particle tracks are reconstructed using the in- 
formation from the DC and the PCI [Iff. The magnetic 
field between the collision vertex and the DC is axial, 
and thus bends particles in the x-y plane. The field is so 
weak at the outer area from the DC that particle tracks 
can be assumed to be straight. A track reconstruction is 
performed in the DC first, and then reconstructed tracks 
are associated with hits in the PCI. The momentum res- 
olution is given by a p /p (%) = 1.3 • p (GeV/c) © 1.0 for 
pions. 

C. Trigger 

The PHENIX experiment has various trigger configu- 
rations to efficiently select many type of interesting rare 
events. This measurement required the coincidence of 
two triggers; a minimum bias (MB) trigger issued by the 
BBC, and a high-energy photon trigger issued by the 
EMCal. 

The MB trigger in p + p collisions requires one charged 
particle in both the north and south sides of the BBC. 
The reconstructed z-vertex is required to be within ± ~ 
30 cm. The efficiency, /mb, of the MB trigger for high-p^ 
QCD scatterings such as jet production is 0.784 ± 0.020, 
which has been determined with the ratio of ir° yields 
with and without the MB trigger requirement. 

The high-energy photon trigger is fired when the sum 
of energy deposits in 4 x 4 EMCal modules (A</> ~ A77 ~ 
0.04) is above a threshold, ~ 1.4 GeV, which varies by 
~ 0.2 GeV area-by-area due to the variations of gain 
and threshold between EMCal modules. Each 4x4 area 
overlaps with others, and thus even when a photon hits 
the edge of a 4 x 4 area the next overlapped 4x4 area can 
gather all energy of the photon. The trigger efficiency is 
almost flat and close to unity above E ~ 2 GeV except 
masked areas due to noise in the trigger electronics. 

III. ANALYSIS METHODS 
A. Outline 

This analysis used 2.3 pb" 1 of data that were taken 
with the MB + high-energy-photon trigger in 2005. In 
addition, ^0.3 pb _1 of data that were taken with the 
MB trigger alone were used for systematic error stud- 
ies. Photons were detected with the EMCal, and charged 
particles were detected with the DC and PCI. Measured 
particles in each PHENIX Central Arm were clustered 
using a cone method to form a 'reconstructed jet' and its 
transverse momentum (py co ). Because of the finite size 
of the acceptance (|f?| < 0.35), the cone size for the parti- 
cle clustering were set to 0.3 at maximum. This is smaller 
than the typical cone size, 0.7 raising two issues: First, 



a jet in an NLO calculation is usually defined with the 
same cone size and compared with the measured jet, but 
this is optimum when both jet energy and cone size are 
large since the jet spread due to hadronization becomes 
significant with small jet energy and cone size. Second, 
such a small cone is more sensitive to quark jets than 
gluon jets since gluon jets are broader and softer than 
quark jets. Because of the situation described above, the 
theory calculation and the simulation evaluations have 
been organized as follows. 

The cross section and the All of inclusive jet pro- 
duction were calculated as a function of jet transverse 
momentum (p^ LO ) within the framework of a next-to- 
leading-order perturbative QCD (NLO pQCD). This cal- 
culation predicted various All's by assuming various 
AG(x) distributions. 

A simulation with the pythia event generator [ij and 
the GEANT detector simulation package [17] was per- 
formed to understand the effects of the detector response, 
the underlying events and the jet-definition difference 
between the measurement and the theory calculation. 
pythia simulates parton-parton hard scatterings in p+p 
collisions at leading order (LO) in a s with phenomeno- 
logical initial and final-state radiation and hadroniza- 
tion. GEANT simulates the acceptance and response of 
the PHENIX detector. We define a jet at the partonic 
level in pythia. The effect of the detector response and 
the underlying events was evaluated as the statistical re- 
lation between the jets defined in pythia and the recon- 
structed jets. We assume p^ Y — py LO within an uncer- 
tainty that will be explained in a later section, and then 
we obtained the relation between the NLO calculation 
and the measurement. 

To confirm that the simulation reproduces well the real 
data in terms of event structure, namely spatial distribu- 
tion of particles in an event, quantities sensitive to event 
structure were measured. Those include particle mul- 
tiplicity, transverse-momentum density, thrust distribu- 
tion and jet-production rate. A comparison was made 
between the real data and the simulation output. 

We derive the predictions of the measured All by con- 
verting the NLO calculation with the relation between 
Pt and py co . A \ 2 test between the measured and 
predicted All's was performed to determine the most- 
probable AG. 

The definitions and relations of jets in this measure- 
ment are summarized in Tab. H] and Fig. [2] 



B. Particle Clustering with Cone Method 

A jet in one PHENIX Central Arm is constructed with 
photons and charged particles detected with the EMCal, 
the DC and the PCI of the Central Arm. A seed-cone 
algorithm, described below, is used for the cluster finding. 
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TABLE I: Definitions of jets adopted in this measurement. 



Reconstructed jet 



jet in pythia 



Hadronic jet made with measurable 
particles after hadronization 
with a cone size of R — 0.3. 
Partonic jet in pythia 
without cone. 



jet in NLO calcula- 
tion 



Partonic jet in NLO pQCD calc. 
with a cone size of 5 — 1.0. 



NLO calculation 



statistical translation 
with PYTHIA+GEANT 
- simulation 



■ measurement 
simulation 



NLO 


<=> 


PY 

Pt 


O 


reco 

Pt 




' cone size: R 



jet (= hard-scattered 
parton) in PYTHIA 




reconstructed 
.et 

measurable 
particles 



FIG. 2: (color online) Relations between the jets defined in 
this measurement. 



1. Event and particle selections 

To select the energy region where the efficiency of the 
high-px* photon trigger is in the plateau, at least one 
photon with pt > 2.0 GeV/c is required in each event. 
This requirement causes a bias towards jets that include 
mostly high-px* ir° , rj, etc. or radiated photons. 

To collect photons from all EMCal hits, a px cut, a 
charged track veto, and an EMCal shower shape cut were 
applied. The pt cut required the pt of each EMCal hit 
to be > 0.4 GeV/c in order to eliminate hits likely to be 
dominated by electronics noise in the detector. It also 
eliminates charged hadron hits because the measured en- 
ergy of minimum ionization particles by PbSc peaks at 
0.25 GeV. and that of tt^ with momentum of 1 GeV/c in 
the PbGl result in a distribution peaked around 0.4 GeV, 
with a broad tail to lower energy. The charged track 
veto reduces charged particle contamination by checking 
whether each EMCal hit has a matched charged track 
within 3ct of their position resolutions. The shower shape 
cut reduces hadron contamination by comparing the frac- 
tion of energy deposits in every EMCal module of a hit 
with the fraction predicted by a model of shower shape. 
This cut eliminates half of hadron hits and statistically 
1% of photon hits. These cuts made the contamination 
of charged and neutral hadrons negligible. 

All charged particles detected with the DC and the 
PCI were required to have pr ranging from 0.4 to 4.0 
GeV/c. Below the lower limit, the acceptance is strongly 
distorted due to a large bending angle and thus becomes 
shifted from that of photons. The upper limit eliminates 
fake high-px tracks which originate from low-px particles 



that are produced from a decay or a conversion in the 
magnetic field. Note that this limit causes a bias towards 
jets that include fewer charged particles. 



2. Cluster finding algorithm 

All particles that satisfy the experimental cuts in one 
arm were used as a seed in cluster finding. Starting with 
the momentum direction of a seed particle as a tempo- 
rary cone axis, we calculated the next temporary cone 
axis with particles which are in the cone. The distance 
between the cone axis (r] C ,(j) C ) and the momentum di- 
rection of each particle {-q 1 , cb l ) is defined as 



kC\2 



(3) 



The cone radius R was set to 0.3, which was about a half 
of the rj acceptance of the detector. The next temporary 
cone axis e nex t is calculated as a vector sum of momenta 
of particles in the cone: 



Cncxt 



p 



\p 



p 



i £ cone 



(4) 



This procedure was iterated until the temporary cone 
axis became stable. 

The cluster finding is done with all seed particles, and 
then each seed particle has one cone and some cones 
can be the same or overlapped. The cone which has the 
largest p™ c ° in an event is used in the event. For mea- 
surements of event structure we also define the sum of 
momenta of all particles in one arm: 



P 



E 



Pi 



(5) 



An evaluation of p T r co without seed has been done us- 
ing a part of the statistics in order to check the effect of 
the use of a seed. Every direction in the (77, (f>) space with 
a step of Srj = 5<p = 0.01 within the Central Arm accep- 
tance has been used as an initial cone direction in each 
event. All steps except the choice of the initial cone di- 
rections is the same as the original algorithm. The yield 
of reconstructed jets with the seedless method was larger 
than that with the seed method by ~20% at p^ c ° = 4 
GeV/c, -10% at p^ co = 8 GeV/c and ~5% at p^ co = 12 
GeV/c. This deviation is compensated in the relation be- 
tween Pt°° and p T Y estimated with the simulation, and 
therefore the p^ c ° difference between the two methods of 
cluster finding is smaller than the deviation above. 



C. Simulation Study 

1. Simulation settings 

The pythia version 6.220 was used. Only QCD high- 
Pt processes were generated by setting the process switch 
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( "MSEL" ) to 1 and the lower cutoff of partonic transverse 
momentum ("CKIN(3)") to 1.5 GeV/c. The parameter 
modification reduces the time for event generation and 
does not affect any physics results in the measured pt 
region, as it has been confirmed by comparing p^ co dis- 
tribution etc. to those without the parameter modifica- 
tion. We call a pythia simulation with these conditions 
'pythia default'. Hadron-hadron collisions have a so- 
called 'underlying event', which comes from the breakup 
of the incident nucleons. The pythia simulation repro- 
duces the underlying event with the Multi-Parton Inter- 
action (MPI) mechanism. The CDF experiment at the 
Tevatron showed that the pythia simulation did not re- 
produce the event structure well and modeled a set of 
tuned parameters called 'tune A' [H, Il8j |. Modified or 
important parameters are listed in Tab. ITT] 



TABLE II: Important or modified (Used) parameters in the 
pythia MPI setting. 



Parameter 


Default 


Used 


Note 


MSTP(81) 


1 


1 


MPI master switch. 


MSTP(82) 


1 


4 


double-Gaussian matter distri- 
bution used. 


PARP(82) 


1.9 


2.0 


turn-off pt for MPI at the ref- 
erence energy scale PARP(89) 




U . o 


U . o 


the fraction of the core Gaus- 
sian matter to total hadronic 
matter 


PARP(84) 


0.2 


0.4 


the radius of the core Gaussian 
matter 


PARP(85) 


0.33 


0.9 


the probability that two glu- 
ons are produced in MPI with 
colors connecting to nearest 
neighbors 


PARP(86) 


0.66 


0.95 


the probability that two glu- 
ons are produced in MPI with 
the PARP(85) condition or as 
a closed loop 


PARP(89) 


1000 


1800 


reference energy scale for the 
turn-off pt 


PARP(90) 


0.16 


0.25 


energy dependence of the turn- 
off pt 


PARP(67) 


1.0 


4.0 


hard-scattering scale p 2 multi- 
plied by this sets the maximum 
parton virtuality in initial- 
state radiation 


MSTP(51) 


7 


7 


CTEQ 5L PDF used. 


MSTP(91) 


1 


1 


Gaussian kr used. 


PARP(91) 


1.0 


1.0 


width of &t distribution. 


PARP(93) 


5.0 


5.0 


upper cutoff for fcr dist. 



We call a pythia simulation with the tune-A setting 
'pythia MPF, although it has been adopted as default 
values in the pythia version 6.226 and later. 

We use the output of the 'pythia default' and the 
'pythia MPF simulations to estimate the effect of the 



underlying event on our measurement. 

The PHENIX experiment has developed its own 
GEANT3-based detector simulator, called the Phenix In- 
tegrated Simulation Application. The absolute scale and 
the resolution of the EMCal energy and the tracking mo- 
mentum have been tuned in the simulation using mass 
distributions of ir° (27), 7r ± , and p ± . 



2. Relation between Pt co and p T 

The pythia+GEANT simulation was used to evaluate 
the effect of the detector response and the underlying 
event on the pJj? co measurement. The px of a jet in 
pythia, which is represented by p^ Y in this paper, should 
be defined so that it is comparable with the theoretical 
jet in order to evaluate the relation between the NLO 
calculation and the measurement. The event- by- event 
transition from the jet in pythia (p^ Y ) to the recon- 
structed jet (pf co ) is simulated to obtain the statistical 
relation between them. 

A jet in pythia is defined as a hard-scattered parton 
that has not undergone final-state parton splits, namely 
particle number 7 or 8 in the pythia event list. A simu- 
lated reconstructed jet is associated with one of the two 
partons by minimizing the angle AR = \J A77 2 + A<fi 2 . 
Figure [3] shows the ratio p^ co /p^ at each p^ co bin, and 
Fig. 0] shows the mean value of the ratios as a function of 
p^ c °. The ratio of the pythia MPI output is -80% on 
average and is larger than that of the pythia default out- 
put due to the contribution from the underlying event. 

The relation between reconstructed jets and jets in 
pythia can be characterized by multiple effects. Some 
particles in a jet can leak from the cone because of the 
limited acceptance, the small cone size and the absence of 
a detector for neutral hadrons. Some particles produced 
by the underlying event can be included in the cone and 
contaminate p^ c °, and thus the ratio Pt co /Pt Y can ex ~ 
ceed one. The p^ Y of events that are in a p^ co bin is 
distributed widely due to the finite pt resolution of the 
PHENIX Central Arm. Because a gluon jet is softer and 
broader than quark jet [l9|, [2(j, the high-pr photon re- 
quirement has lower efficiency for gluon jets. Therefore 
the ratio of p^ c ° to p^ Y for gluon jets is smaller than 
quark jets on average. 

Figure [5] shows the relative yields of quark+quark 
(q + q), quark+gluon (q + g) and gluon+gluon (g + g) 
subprocesses as a function of p^ Y at each p^ co bin. Fig- 
ure [5] shows the fraction of g + g, q + g and q + q sub- 
processes as a function of p™ c °. These were evaluated 
with the simulation. As explained above, the gg subpro- 
cess is suppressed in this measurement. The dominant 
subprocess is q + g throughout the p^ co range. 
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FIG. 3: (color online) Distributions of the ratio Pt c °/Pt Y 
evaluated with (dashed black) pythia default and (solid 
green) pythia MPI. 



FIG. 5: (color online) The relative yields of q + q, q + g and 
g + g subprocesses in the pythia+GEANT simulation. The 
results with all the subprocesses combined are also shown. 
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FIG. 4: (color online) The mean value of ratio Pt^/Pt* as 
functions of py co . 



3. Relation between p? Y and Pr h ° 



The cross section and the All of inclusive jet pro- 
duction in |?7| < 0.35 at y/s = 200 GeV were calculated 
within the NLO pQCD framework with the CTEQ6M 
unpolarized PDF under the Small Cone Approximation 
(SCA) [HHH. We adopted a cone size of 8 = 1.0 for 
reasons that will be explained in a later section. Figure 



[7] shows the cross section calculated with three factor- 
ization scales, /i = pt, %Pt and pt/2 in NLO pQCD. 

The p^ LO needs to be connected with p^ Y in order to 
evaluate the relation between the NLO calculation and 
the measurement, where the relation between p^ Y and 
preco wag obtained from the pythia+GEANT simulation. 
We assume p^ Y — p^ LO , and thus the relation between 
the jet in pythia and the measurement can be inter- 
preted as the relation between the NLO calculation and 
the measurement. However the definition of p^ Y and 
p^ LO has a discrepancy, and they become close to each 
other only as the cone half-aperture (S) in the theory 
becomes large. Therefore we set 6 to 1.0, which is the 
upper limit where the SCA is applicable, and evaluated 
the discrepancy between p^ Y and p^ LO with S = 1.0 as 
described later. Moreover, the cone size of the jet in the 
NLO calculation needs to be larger than the acceptance 
of the PHENIX Central Arm so that one jet per cen- 
tral arm per event can be reconstructed and connected 
with the jet in the NLO calculation. This has been also 
satisfied with the use of S — 1.0. 

Note that the cone size in theory and measurement are 
different parameters and the difference is compensated 
for with the pythia simulation; the former is related to 
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FIG. 6: (color online) Subprocess fractions of reconstructed 
jets as functions of pr- It was evaluated with the pythia MPI 
and GEANT simulation. It should be noted that the gluon- 
quark reaction is the dominant reaction in all the momentum 
region from 4 to 12 GeV/c. 
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FIG. 7: Unpolarized jet cross section at a pseudorapidity 
\r)\ < 0.35 with a cone half-aperture 5 = 1. It was calculated 
at NLO under the SCA with three factorization scales, /i = pr 
(solid line), 2px (lower dashed line) and Pt/2 (upper 
dashed line). 



the angle between two splitting partons and the latter is 
related to the angle between stable particles. 



that the scales of p™ conc and p^ LO are the same. Then 
the difference between conG and p T Y , which can be 
evaluated using pythia, is considered to be the differ- 
ence between p^ LO and p T Y ■ 

Figure [S] shows distributions of the fraction 
pin cone^PY three typical p T Y bins. This indi- 
cates that the px scales of the two jet definitions have 
a 10% difference on average in the pt range of these 
measurements. Therefore the uncertainty due to the 
jet-definition difference between pythia and the NLO 
calculation with S = 1.0 has been assigned 10% in px 
scale. 
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FIG. 8: (color online) Distributions of the fraction 
pin concypPY evama t ec ! w ith a pythia simulation at three typ- 
ical py Y bins. 



5. Reproducibility Check 

Figure [9] shows the distribution of p^ c ° measured with 
the clustering method described above. The simulation 
outputs have been normalized so that they match the real 
data at p T ~ 8 GeV/c. The slope of the pythia MPI out- 
put agrees better with that of the real data, where that 
of the pythia default output is less steep. The relative 
yield between the real data and the PYTHIA MPI output 
is consistent within ±10% over five orders of magnitude. 



4- Uncertainty due to difference in jet definitions 

The uncertainty due to the jet-definition difference be- 
tween the pythia and NLO calculations with 5=1.0 has 
been evaluated using the difference between two jet def- 
initions in pythia. One definition is the jet in pythia 
defined above. The other assumes a cluster of partons 
with a cone size of S = 1.0 in pythia, where partons 
originating from the underlying event are excluded. For 
the latter definition the jet px is denoted p™ cono . Since 
pin cone anc j pNLO are (Fefined similarly, i.e. both at the 
partonic level and with the same cone size (5, we assume 



Figure [10] shows distributions of the fraction 
ptrigri J pveco ^ wnere p^ lg ^ ig p T f the trigger photon. The 
lower cutoff of the distributions is due to the minimum 
Pt of the trigger photon (> 2 GeV/c). The rightmost bin 
(Py 1S7 /Px c ° ~ 1) contains events in which only a trigger 
photon exists. Such events can occur by the limited ac- 
ceptance, by the EMCal masked area (particles except a 
trigger photon in jet are not detected) , by EMCal noise or 
by direct photon events. The difference between the real 
data and the simulation outputs in the rightmost bin may 
indicate that these effects are not completely reproduced 
by the simulation, but the difference is small (< 5%) and 
negligible in comparison with other uncertainties. 
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FIG. 9: (color online) Reconstructed-jet yields as a function 
of Pt co . The red, black and green points correspond to the 
real data, the pythia default output and the pythia MPI 
output, respectively. The simulation outputs have been nor- 
malized so that they match the real data at pr = 8 GeV/c. 
The ratio of the yields between the simulations and the real 
data is shown at bottom. 
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FIG. 10: (color online) The fraction of pt of the trigger 
photon in each p™ co . 



IV. RESULTS AND DISCUSSIONS 
A. Event structure 

1. Multiplicity 

Multiplicity is defined as the number of particles which 
satisfy the experimental cuts in one event. Figure [TTT a) 
and (b) show the mean value of multiplicity in the Central 
Arm vs p^ m and in the cluster vs p^ c ° . The multiplici- 
ties in the arm and in the cluster of the simulation out- 
puts agree, on the whole, with that of the real data. The 
pythia MPI output is larger than the pythia default 
output as expected, and the real data are closer to the 
pythia default output. On the other hand, the p^ c ° dis- 
tributions (Fig. [5]) shows better agreement between the 
real data and the PYTHIA MPI output. This indicates 
that the pythia MPI reproduces the sum of px of parti- 
cles well, which is less sensitive to particle fragmentation 
process, while it does not reproduce the particle multi- 
plicity very well. The reproducibility of the summed pt 
is checked in measurements described later. 

Figure fTTTc) and (d) show the ratio of charged-particle 
multiplicity to photon multiplicity in the Central Arm 



and in the cluster. The real data lies below the pythia 
default and MPI results for both multiplicities. This in- 
dicates that the effect of the underlying event in the ra- 
tios is small, and the difference between the real data and 
the PYTHIA results is mainly caused by the imbalance be- 
tween photons and charged particles in jet. Figure ITTT e) 
and (f ) show the ratio of the sum of charged-particle pt 
to the sum of photon px ■ These have the same tendency 
as the multiplicity ratios. 



2. Transverse momentum density 
The pt density, T)p T (A<p), is defined as 

V PT (A^^/j- PTi) , (6) 

\ 9 im[A& / cvont 

where Acj) is <f> angle with respect to the direction of a 
trigger photon in event, 5(f) is an area width in </> direc- 
tion, and pt% is transverse momentum of i-th particle in 
event. The px density means the area- normalized total 
transverse momentum in an area of 5<f> x Srj at a distance 
Acj) from trigger photon, where 8rj is the width of the 
Central Arm acceptance. 
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FIG. 11: (color online) (a): Mean multiplicity in the Central 
Arm vs p^ m . (b): Mean multiplicity in the cluster vs p r f co . 
(c): The ratio of charged multiplicity to photon multiplicity 
in the Central Arm. (d): Same as (c) but in the cluster, (e): 
The ratio of charged pr to photon pr in the Central Arm. 
(f): Same as (e) but in the cluster. 



We name the region at A0 < 0.7 rad the 'toward' re- 
gion and the region at A<fi > 0.7 rad the 'transverse' 
region. Since particles from a jet are concentrated along 
the jet direction, the T> PT in the transverse region is sen- 
sitive to the underlying event. 
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FIG. 13: (color online) pr density, T)p T = dSiPTi/d(fi 
(GeV/c/rad), in each pj?™ 1 bin. Trigger photons are included 
in the leftmost points. 

Figure [13] shows the V PT distributions for each p?^ m 
range. In the "toward" region, the simulation outputs 
agree well with the real data. It shows that the shape 
of jets produced by the simulation is consistent with the 
real data. In the "transverse" region, the PYTHIA default 
output is generally smaller than the real data. This is 
an indication that the PYTHIA default does not contain 
sufficient total pt of soft particles from the underlying 
event. The PYTHIA MPI output agrees with the real data 
well. 



FIG. 12: (color online) Measurement condition of the pr 
density. The arc and the x mark represent one Central Arm 
and the collision point in the beam view. 

As illustrated in Fig. [12l to avoid the effect of the 
PHENIX Central Arm acceptance in the calculation of 
2? PT , we limited the <fi direction of the trigger photons 
to less than 20° from one edge of the PHENIX Central 
Arms, and we did not use photons and charged particles 
which were in the <p area between the trigger photon and 
the near edge. With this method the T> PT distribution 
is not affected by the finite acceptance of the PHENIX 



3. Thrust distribution in PHENIX Central Arm 

We evaluated the thrust variable defined in the CERN- 
ISR era with particles in one PHENIX Central Arm 
(At/ = 0.7, A(f> = 90°): 



Tph = max 



E t \Pi ■ u \ _ Ei \Pi ■ gj 



p 



Ei \Pi 



E< \Pi 



(7) 



(8) 
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where u is a unit vector which is called the thrust axis 
and is directed to maximize T, and p i is a momentum 
of each particle in one arm. If only particles in a half 
sphere in an event are used, Tpu can be written as the 
right-side formula in Eq. [7] 

The distribution of Tpn of isotropic events in the 
PHENIX Central Arm acceptance for each p™ m bin was 
simulated with the following method. First, the cross 
section of inclusive particle production is assumed to 
be proportional to exp(— 6 pxiGeY /c)) and is indepen- 
dent of rj and <f>. Second, the same cuts as the experi- 
mental conditions are applied numerically: the geomet- 
rical acceptance (\n\ < 0.35, A<fi = 90°), the momen- 
tum limit (j>t > 0.4GeV/c), and one high-pr particle 
{p T > 2.0GeV/c). Third, the distribution of T PH of 
isotropic events was calculated for each number of par- 
ticles in one event (f n (T) for n = 1, 2, 3, ... ). The Tpn 
distribution of n = 2 events is particularly steep. Thus 
we applied a cut of n > 3 in the Tph measurement. The 
fr is evaluated as the sum of / n (T)'s weighted by the 
probability (e„) that the number of particles per event is 



/(T) = ^e„/„(T) 



(9) 



where e„ was derived from the real data. 

Fi gure 1141 shows the Tph distribution in each p?^ nl 
range. The pythia MPI output agrees with the real data 
well. The pythia default has a steeper slope, which indi- 
cates that the number of particles in the vicinity of jets 
in the PYTHIA default is insufficient. In the real data, 
the pythia default output and the pythia MPI output, 
the Tpn distribution becomes sharper as p^ m increases. 
This is due to the fact that the transverse momentum 
(jt) of a jet is independent of its longitudinal momen- 
tum and is almost constant. 

If the real data includes a contribution from non-jet 
(isotropic) events, the Tph distribution of the real data 
is a mixture of the distribution of the simulation output 
and the distribution of the isotropic case. The contri- 
bution from non-jet events can be judged to be negligi- 
ble because the pythia MPI output reproduces the data 
even though it docs not have isotropic events. 
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FIG. 14: (color online) Tph distribution in each py um bin. All 
distributions have been normalized so that their areas were 
equal to one another. The purple lines are the distributions 
of isotropic events in the acceptance of the PHENIX Central 
Arms, which are evaluated with Eq. [9] 



the high-pr photon trigger, respectively; N^. eco is the 
reconstructed-jet yield in a z-th p r rp c ° bin. The high- 
Pt photon trigger efficiency f p h was estimated to be 
0.92 ±0.02, where the inefficiency is caused by the 10% of 
the EMCal acceptance where the trigger was disabled due 
to electronics noise. The inefficiency is slightly smaller 
than the disabled acceptance because a particle cluter 
can contain multiple high-p^ photons. 





B. Jet production rate 

1. Evaluation method (measurement) 

The jet production rate y, namely the yield of recon- 
structed jets per unit luminosity, is defined with mea- 
sured quantities as 

N* 

y ee - — , (io) 

L ■ /MB ' Jph 

where L is the integrated luminosity; /mb and f p h are 
the efficiencies of the MB trigger (see Sec. Ill CI) and 



2. Evaluation method (prediction) 

On the other hand, the variable y is expressed with 
theoretical and simulation quantities as 

y = f J ' e trig+acc ' 3 / 't?ieo' (H) 

3 

where the label i and j are the indices of p^ co and p^, LO 

bins, respectively. The y 3 t heq ' IS a production rate 
within |rj| < 0.35 in a j-th p^ LO bin, which is theoreti- 
cally calculated. The £ , trig+acc is a high-p^-photon trigger 
efficiency and acceptance correction, which is evaluated 
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with the PYTHIA+GEANT simulation. The ^^^^{^0 
is a yield of jets that include a high-p^ photon within 
|ry| < 0.35. The / IJ is the probability that a jet within 
a j-th p^ LO bin is detected as a reconstructed jet within 
a i-th p!jf co bin. This method uses the relative p!j? co dis- 
tribution in each bin and thus the slope of the p^ Y 
distribution in the simulation does not affect the result 

The correction factor e J trig+acc is a fraction, whose nu- 
merator is the number of events in which at least one 
photon with p T > 2 GeV/c is detected, and whose de- 
nominator is the number of events in which jets are in 
\rj\ < 0.35. The condition ll p T > 2 GeV/c" in the numer- 
ator corrects a high-pr photon efficency, i.e. the probabil- 
ity that a high-pT photon in jets must be detected with 
the EMCal. The condition "\rj\ < 0.35" in the denomina- 
tor and the absence of it in the numerator corrects an ac- 
ceptance for jets, i.e. the fact that a part of reconstructcd- 
jets does originate from jets with \r/\ > 0.35. Figure [T5l 
shows 



trig+acc 



pythia default and MPI simulations. 
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FIG. 15: (color online) The correction factor e 3 trig+acc for 
high-pr photon trigger efficiency and acceptance effect. The 
pythia default (black) and the pythia MPI setting (green) 
were used. 



To estimate a systematic error related to the simula- 
tion reproducibility of high-p^ photon, we evaluated, in 
both the real data and the simulations, the ratio (r) of 
the reconstructed-jet yields in the high-pr photon trig- 
gered sample to that in the MB triggered sample. The r 
of the pythia MPI output is 5% at p^ co = 4 GeV/c and 
50% at p!j5 co = 12 GeV/c, and is consistent with that of 
the real data within ±10%. Therefore a 10% error was 
assigned to the jet production rate calculated with the 
PYTHIA MPI simulation. The r of the pythia default 
output is smaller by 20-30% than that of the real data. 



3. Result 

Figure [T5] shows the jet production rate. The main 
systematic errors are listed in Tab. IIIII The main uncer- 



tainties of the measurement are the BBC cross section 
and the EMCal energy scale. These errors are fully cor- 
related bin-to-bin. The error on the EMCal energy scale 
includes both the change of px of individual photons and 
the change of the threshold of the high-p^ photon require- 
ment. In comparing the measurement and the calcula- 
tion, the 10% pt scale uncertainty of the jet definitions 
in the pythia simulation and the NLO pQCD theory 
makes a 30% error at low pt or 70% at high py, and is 
the largest source. The uncertainty of the renormaliza- 
tion and factorization scales in the NLO jet production 
cross section makes a 30% error. The calculation with 
pythia MPI agrees with the measurement within errors 
over the measured range 4 < p!j? co < 15 GeV/c. 



as a function of py LO estimated with the g i» ! 
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FIG. 16: (color online) (top) Reconstructed jet yield and 
(bottom) the ratio of the real data to the calculations, (red 
points) Real data with the (gray band) experimental system- 
atic error, (black curves) pythia MPI calculation with theory 
factorization scales of (solid curve) pr, (upper dashed curve) 
Pt/2, and (lower dashed curve) 2pr- (dotted curves) Varia- 
tion caused by 10% pr scale uncertainty around the calcula- 
tion. Statistical uncertainties are smaller than the size of the 
points. 

The result with pythia default is smaller than the 
result with pythia MPI by 50% at p^ co = 4 GeV/c, 
by 35% at p^ co = 9 GeV/c and by 20% at p^ co = 14 
GeV/c. It can be fully explained by the difference visi- 
ble in Fig. H"5l between pythia default and PYTHIA MPI. 
According to the comparisons of the event structure, 
pythia MPI reproduces the spatial distribution of parti- 
cle momenta in one event much better than the pythia 
default. Therefore, for the jet production rate evaluated 
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TABLE III: Main systematic errors of the jet production rate. 



Source 


Size 


Size on rate 


Measurement 






Luminosity 


9.7% 


9.7% 


EMCal energy scale 


1.5% 


7-6% 


Tracking momentum scale 


1.5% 


0-3% 



Calculation 



Jet definition 10% in p T 30-70% 

Jet shape & underlying event 50-20% 

High-pT photon fragmentation 10% 

Simulation statistics 2-5% 



and the reconstructed z-vertex is within 30 cm of the 
collision point. 

The beam polarizations were measured with the pC 
and H-jet polarimeters [13, HH at the 12 oclock inter- 
action point on the RHIC ring. One of the colliding 
beam rotating clockwise is called "blue beam", and the 
other rotating counterclockwise "yellow beam". The 
luminosity-weighted-average polarizations are 50.3% for 
the blue beam and 48.5% for the yellow beam. The sum 
of statistical and systematic errors on (Pb){Py) is 9.4%. 

2. Evaluation method (prediction) 



with PYTHIA MPI simulation, the error due to possible 
insufficient tunings of pythia MPI should be smaller 
than the difference of the jet production rate between 
the PYTHIA MPI simulation and the PYTHIA default sim- 
ulation. 



C. Double helicity asymmetry All 



Evaluation method (measurement) 



All is expressed with measured quantities as 



A LL = 



1 (N++ + N-) - R(N+- + N-+) 
\P b \\Py\ (N++ + V__) + R(N + _ + N_+) 



R = 



L+++L- 



L. 



(12) 



(13) 



where N++ etc. are reconstructed-jet yields with collid- 
ing proton beams having the same (H — h or ) and op- 
posite (H — or — h) helicity; Pb and Py are the beam 
polarizations; R is the relative luminosity, i.e. the ratio 

of the luminosity with the same helicity (L ++ + L ) 

to that with the opposite helicity (L+- + L |_). All is 

measured fill-by-fill and the results are fit to a constant, 
because the beam polarization and the relative luminos- 
ity are evaluated fill- by-fill to decrease systematic errors. 
The average fill length was about five hours. The inte- 
grated luminosity used was 2.1 pb^ 1 . It is 0.1 pb -1 less 
than the statistics used in the production rate measure- 
ment because the data with bad conditions on the beam 
polarization were discarded. 

The relative luminosity at PHENIX was evaluated 
with the MB trigger counts (A^g and -V^g) as R = 
^mb/^mb- A possible spin dependence of MB-triggered 
data causes an uncertainty on the relative luminosity. 
The error has been checked by comparing the relative 
luminosity with another relative luminosity defined with 
the ZDCLL1 trigger counts. The ZDCLL1 trigger is fired 
when both the north ZDC and the south ZDC have a hit 



Polarized/ unpolarized cross sections of jet production 
for every subprocess (q + q, q + g and g + g) were cal- 
culated at NLO based on the SCA with a cone size of 
8 = 1.0. The polarized cross sections were calculated us- 
ing various AG(x) in order to compare the measured All 
with various predicted All 's and find the most-probable 
AG(x). Figure [TTl shows the distributions of the AG(x) 
used, and the integrated values are 



dxAG(x, n 2 -- 

' -1.24 (AG 
-1.05, 
-0.90, 
-0.75, 
-0.60, 
-0.45, 
-0.30, 
-0.15, 



0.4GeV 2 

:-G), 



'0 (AG = 0), 

0.24 (GRSV-std), 

0.30, 

0.45, 

0.60, 

0.70, 
[1.24 (AG = G) 



(14) 



Each AG(x) (except the GRSV-std, the AG = G input, 
the AG = input and the AG = — G input) have been 
obtained by refitting the GRSV parameters to the DIS 
data which were used in the original GRSV analysis [26[ . 
It is noted that the DIS data used in GRSV are the data 
up to the year 2000 and thus are much less than that 
used in the updated analysis, DSSV [27J, for example. 
The polarized PDF in the GRSV parameterization is of 
the form: 

Af(x,vo 2 ) = N f x a ?(l - xfff(x,ti 2 )GRV, (15) 

where / is u, d, q or G; /io 2 = 0.4 GeV 2 is the initial 
scale at which the functional forms are defined as above; 
f(x, Unf)GRV is the unpolarized PDF of the GRV98 anal- 
ysis [28j; Nf, af and j3f are free parameters. In the refit 
of the DIS data, the integral value of AG(x) from x = 
to 1 was fixed to its particular value listed above, and the 
shape of AG(x) and the quark-related parameters were 
made free. The \ 2 of the refitting to the DIS data is 170 
for the 209 data points [26| when the integral of AG is 
at the initial /x 2 , for example. In the remainder of this 
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FIG. 17: (color online) Assumed gluon distribution functions 
at jj 2 = 1 GeV 2 . The integral J Q dxAG(x) of each distribution 
at the initial scale /i 2 = 0.4 GeV 2 is, from bottom to top at 
x = 0.15, -1.24 (AG = — G), -1.05, -0.90, -0.75, -0.60, -0.45, 
-0.30, -0.15, (AG = 0), 0.24 (GRSV-std), 0.30, 0.45, 0.60, 
0.70 and 1.24 (AG = G). 



paper we concentrate on investigating the y 2 of the six 
data points of the reconstructed-jet All- 

The various AG(x) above were evolved up to a scale 
(i of every event in the All calculation. The All of 
every subprocess (A q L + L q , A q L + L 9 and A 9 L + L 9 ) can be derived 
as functions of p^, from the unpolarized and polarized 
cross sections. The pythia+GEANT simulation produces 
the relative yields of every subprocess (n 9+9 (py ,Py co ), 

n g+5(^NLO )p reco-) and n g+g ( p NLO ^ p rcco^ ag ghown in 

Fig. EJ A"£ c L °{p^ co ) is calculated as a mean of A q L +q , A q L + L 9 
and A{ + L 9 weighted by the fractions of events: 



a rcco / rcco 
A LL (Pt 



I 



isub 



i <8ub (p^ LO ,p5? co ) • At L h (PT LO ) 



zsub 



(16) 



,Pt 



where isub is q + q, q + g and g + g. As an estimation of 
systematic errors, the slope of jet yields and the fraction 
of subprocesses were compared between the theory cal- 
culation and the pythia simulation. Note that both the 
slope and the fraction that we compared have not been 
biased by the high-p^ photon and the small cone, since 
the theory calculation cannot provide biased values. The 
variations of A^£ caused by both the slope difference 
and the fraction difference are negligible in comparison 
with other errors. 



Result 



Figure [TS] shows measured and four predic- 

tion curves. Table IIVI shows the values of measured 
^jroco_ rp^g measu red All is consistent with zero, as the 



X 2 /n.d.f. between the data points and zero asymmetry 
(All = 0) is 1.3/6. The systematic error of the relative 
luminosity is much smaller than the statistical error on 
All and is negligible. On the prediction curves the sys- 
tematic error related to the fractions of subprocesses are 
smaller than the 10% pr scale uncertainty by roughly an 
order of magnitude. Therefore it is not included in this 
plot. 
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FIG. 18: (color online) Reconstructed-jet All as a function 
of pf c °- (red points) Measurement with statistical error bars, 
(black lines) Calculation based on four AG(x) functions and 
the pythia MPI + GEANT simulation. 



TABLE IV: Measured reconstructed-jet All- 



Py co range and mean (GeV/c) 


All 


stat error 


4-5, 4.42 


-0.0014 


0.0037 


5-6, 5.43 


-0.0005 


0.0059 


6-7, 6.43 


0.0058 


0.0089 


7-8, 7.44 


0.0034 


0.0132 


8-10, 8.79 


0.0077 


0.0152 


10-12, 10.81 


-0.0181 


0.0282 



It has been confirmed with a "bunch shuffling" method 
that the size of the statistical errors assigned is appropri- 
ate. In this method, the helicity of every beam bunch was 
newly assigned at random and A^ c £° was evaluated again. 
Repeating this random assignment produced a large set 
of A T £^° values. Its mean value should be of course zero 
and was confirmed in this exercise. Its standard devia- 
tion indicates the size of the statistical fluctuation, and 
was consistent with the statistical errors assigned. The 
point-to-point variance seems smaller than the statisti- 
cal errors of the data points, but we could not find any 
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unrecognized cause such as a statistical correlation. We 
conclude that the small variance of the data points hap- 
pened statistically despite its small probability. 

As a systematic error check, the single spin asymmetry 
Al was measured. It is defined as 



At EE 



1 N + — R N_ 
PN+ + RN- 



R = 



(17) 



where N + and iV_ are reconstructed-jet yields with one 
colliding proton beam having the positive and negative 
helicity, respectively; P is the beam polarization; R is 
the relative luminosity, i.e. the ratio of the luminosity 
with the positive helicity (£+) to that with the negative 
helicity As the jets are produced via the strong 

force, Al must be zero under the parity symmetry. Thus 
any non-zero value indicates systematic errors. 

Figure shows measured A^. Al was measured for 
the polarization of one colliding beam while the other 
beam was assumed to be unpolarized. No significant 
asymmetry was observed. 
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FIG. 19: (color online) Jet Al as a function of p T f c ° ■ The 
blue and green points are the results using the polarizations 
of the blue beam and the yellow beam, respectively. The red 
points are the averages of the blue and green points. 



4- Constraint on AG 

To determine the range of x g i uon probed by this mea- 
surement, the PYTHIA MPI simulation without GEANT 
was used to obtain event- by- event x g i uon (one value per 
q-g scattering event, two values per g-g, or none per q-q) 
and also n 2 . Figure [501 an d HH show the distributions of 



Xglv 



and fj. , respectively. The x g i uon value where the 



yield is half maximum is 0.02 at the lower side of the 
"4 < p™ co < 5" distribution and 0.3 at the upper side of 
the "10 < p5? co < 12" distribution. Therefore we adopt 
a range of 0.02 < x g i uon < 0.3 as the range probed by 
this measurement. Table IVl shows the integral of AG(x) 
at the measured x g i uon range, below the range and above 
the range. The measured x g i uon range includes ~70% of 



distributions in all the four GRSV models shown. With 
the same procedure, the /i 2 range probed was estimated 
to be 5 < a* 2 < 300 GeV 2 . 



TABLE V: Partial integral of AG(x) at /x 2 = 1 GeV 2 



Model 


J dxAG(x) at each x ranj 


;e 




lCT 4 -0.02 


0.02-0.3 


0.3-1 


10" 4 -1 


AG = -G input 


-0.406 


-1.09 


-0.208 


-1.71 




(24%) 


(64%) 


(12%) 




AG = input 


0.00808 


0.0644 


0.00869 


0.0812 




(10%) 


(79%) 


(11%) 




GRSV-std 


0.0684 


0.258 


0.102 


0.427 




(16%) 


(60%) 


(24%) 




AG = G input 


0.427 


1.22 


0.226 


1.87 




(23%) 


(65%) 


(12%) 
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FIG. 20: (color online) Distributions of x g i uon in events that 
include a reconstructed jet with 4 < p r f co < 12 GeV/c. 
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FIG. 21: (color online) Distributions of fj, 2 in events that 
include a reconstructed jet with 4 < py co < 12 GeV/c. 

Figure [2U shows the y 2 between the 6 data points 
and the prediction curves as a function of the integral 
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Jo^ dxAG(x, n 2 = 1) for each prediction curve. The 
value of fj, 2 (= 1 GeV 2 ) has been arbitrarily chosen in 
order to show the value of the AG integral in horizon- 
tal axis. Actual /j, 2 used in the All calculation varies 
depending on jet pr- 
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FIG. 22: (color online) \ 2 between the measured All and 
the calculated All as a function of the integrated value of 
AG(x). 

The minimum of the \ 2 is ~ 1-5 at AG = 0.07, namely 
the GRSV AG = input. The 95% and 99% confidence 
limits are where the x increases from the minimum by 
4 and 9, respectively. We obtained 



/■0.3 

1.1 < / AG 

JO. 02 



GRSV 



(x, fx 2 = 1) < OA (18) 



at 95% confidence level and 



The MPI-enhanced pythia simulation agrees well with 
the real data in terms of the event structure: the multi- 
plicity of photons and charged particles, the pt density 
as a function of the azimuthal angle from trigger photon, 
and the thrust in the PHENIX Central Arm. A small 
difference in the intra-jet structure, namely the fractions 
of photons and charged particles in jets, was observed 
as shown in Fig. UTT c) to (f). Nevertheless, the simula- 
tion well reproduces the shape of jets and the underlying 
event at this collision energy. 

In the measurement of jet All, measured particles 
were clustered by the seed-cone algorithm with a cone 
radius R = 0.3. The relation between p^ LO and p^ co was 
evaluated with pythia and GEANT. The jet production 
rate was measured and satisfactorily reproduced by the 
calculation based on the NLO pQCD jet production cross 
section and the simulation. The jet All was measured 
at 4 < p r rf co < 12 GeV/c. The main systematic errors are 
a Pt scale uncertainty of 10% and a beam polarization 
uncertainty of 9.4%. The x g i uon range probed by this 
jet measurement with 4 < p^ c ° < 12 GeV/c is mainly 
0.02 < x < 0.3 according to the simulation. The mea- 
sured All was compared with the predicted values based 
on the GRSV parameterization, and the comparison im- 
posed the limit -1.1 < f ° 3 2 dxAG GRSV (x, y 2 = 1) < 0.4 

at 95% confidence level or f°* 2 dxAG GRSV (x, fi 2 = 1) < 
0.5 at 99% confidence level. The theoretical uncertainties 
such as the parameterization of the polarized PDFs were 
not included in this evaluation. 
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0.3 



AG 



GRSV 



(x,n 2 



1) < 0.5 



(19) 



0.02 



at 99% confidence level. In the assumptions of the 
present approach, the error correlations between the 
normalization parameter and the shape parameters in 
AG(x) are not included. Also the fact that the shape 
of the polarized PDFs is parameterized into Eq. [15] may 
cause additional uncertainty in AG(x). 



V. CONCLUSION 

We measured the event structure and the double hc- 
licity asymmetry (All) in jet production at midrapidity 
(|ry| < 0.35) in longitudinally polarized p + p collisions 
at yfs = 200 GeV were measured. The main motiva- 
tion is to use this complementary approach to inclusive 
measurements to better understand the contribution of 
the gluon spin (AG) to the proton spin. Because this 
measurement of All observes a larger fraction of the jet 
momentum, it reaches higher px and thus higher gluon 
x. 
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